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ABSTRACT 

The  purpose  of.  this  study  is  to  investigate  the 
influence  of  non-environmental  and  small  scale  environ- 
mental acoustic  factors  on  the  variation  of  transmission 
loss  within  areas  of  the  ocean  which  exhibit  homogeniety  of 
large  scale  environmental  factors.   The  method  employed  was 
to  identify  the  influencing  and  non-influencing  factors, 
pool  over  the  non-influencing  factors  if  they  existed,  and 
investigate  the  nature  and  magnitude  of  the  effect  of  the 
influencing  factors  on  the  variation  of  transmission  loss. 

The  results  generally  indicated  that  source/receiver 
depth,  the  mixed  layer,  geographic  location,  and  transmit- 
ting range  and  frequency  significantly  influenced  the 
variation  of  transmission  loss. 
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I.   INTRODUCTION 

In  recent  years,  with  the  advent  of  more  sophisticated 
SONAR  systems,  there  has  been  an  increasing  emphasis  placed 
on  the  study  of  the  acoustic  properties  of  the  ocean.   To 
be  able  to  predict  the  capabilities  of  a  particular  SONAR 
system,  whether  it  be  passive  or  active  in  nature,  requires 
a  complete  understanding  of  the  effect  of  the  local  acous- 
tic environment  on  the  transmission  of  sound  through  the 
water.   Both  theoretical  and  observational  techniques  have 
been  utilized  in  an  effort  to  determine  the  environmental 
and  non-environmental  factors  which  effect  the  transmission 
of  sound  and  then  to  characterize  this  transmission  as  a 
function  of  those  factors. 

One  study  of  this  nature  was  the  AMOS  project  performed 
under  the  sponsorship  of  the  U.  S.  Navy  in  the  early  1950' s 
(ref.  5).   In  this  study,  the  behavior  of  sound  transmission 
was  measured  along  with  local  environmental  characteristics 
over  a  wide  range  of  environmental  and  non-environmental 
factors.   On  the  basis  of  these  observations,  transmission 
loss  prediction  equations  were  developed  as  a  function  of 
certain  factors  such  as  source  and  receiver  depths,  range  of 
transmission,  mixed  layer  depth,  etc.   However,  these  equa- 
tions are  quite  insensitive  to  small  scale  environmental 
characteristics  such  as  inhomogenieties  in  the  thermal 
structure  and  irregularities  at  the  air-sea  and  sea-earth 


interfaces.   As  a  consequence,  transmission  losses  observed 
operationally  are  often  quite  variable  with  respect  to  the 
predicted  losses  on  the  basis  of  these  equations. 

This  study  is  an  extension  of  the  work  performed  in  the 
AMOS  experiment  in  that  it  investigates  the  nature  of  the 
small  scale  environmental  factors  which  influence  the 
observed  temporal  and  spacial  variability  in  acoustic  pro- 
pogation.   The  characterization  of  this  variability  as  a 
function  of  these  factors  will  be  useful  in  a  variety  of 
applications.   Some  examples  are: 

1)  Establishing  some  design  specifications 

for  new  SONAR  systems. 

2)  Establishing  SONAR  detection  probabilities 

for  shipboard  use. 

3)  Simulation  of  environmental  conditions  for 

use  with  detection  simulation  models. 
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II.   A  MODEL  OF  THE  ACOUSTIC  ENVIRONMENT 

Within  sufficiently  small  areas  of  the  ocean, "the  large 
scale  environmental  factors  (i.e.,  thermal,  density  and 
salinity  profiles)  may  be  considered  spacially  consistent. 
Such  areas  may  also  be  considered  temporally  consistent  for 
short  durations  of  time.   Since  general  acoustic  behavior 
is  a  function  of  these  large  scale  environmental  conditions, 
such  an  area  may  be  considered  acoustically  homogeneous. 
This,  however,  does  not  imply  that  the  acoustic  behavior  is 
independent  of  small  scale  temporal  and  spacial  inhomogeni- 
eties  within  this  large  scale  environmental  structure. 
Hence,  the  environment  may  be  characterized  as  a  generally 
homogeneous  medium  which  dictates  the  large  scale  acoustic 
behavior,  but  one  having  small  temporal  and  spacial  inhomo- 
genieties  which  cause  variations  in  acoustic  behavior. 

On  the  basis  of  this  characterization,  sound  transmission 
loss  may  be  modeled  as  a  stationary  stochastic  process 
whose  mean  is  a  function  of  the  large  scale  environmental 
factors,  and  whose  autocovariance  function  is  a  function  of 
the  small  scale  environmental  factors.   Denote  this  process, 
where   X(t)   describes  the  transmission  loss,  by   X(t);t>0. 
It  is  convenient  to  characterize  this  process  by  its  mean 

/c  =  E(X(t)) 

and  by  its  autocovariance  function 


C(  A  t)  =  E((X(t)  -  yU)  (X(t  +  A  t)  -/<.)). 

It  is  then  possible  to  study  the  nature  of  the  variability 
which  occurs  in  transmission  loss  due  to  small  scale 
environmental  aspects  by  considering  the  environmental  and 
non-environmental  factors  which  effect  C(  A  t) . 

Two  particular  parameters  which  characterize  this 
variability  are   C(0),  the  variance  of  the  process  and  hence 
the  amount  of  variability  that  would  be  experienced  within 
an  acoustically  homogeneous  region  for  a  given  source  and 
receiver  depths,  and   At  ,  the  minimum  value  of   At   for 
which   C(At)  =0,  which  represents  the  shortest  time  inter- 
val between  independent  observations.   These  two  parameters, 
C  (O)   and   At*,  may  be  determined  for  a  given  process  if 
that  process  is  observed  either  continually  or  at  discrete 
increments  of  time.   Then  by  observing  such  processes  under 
varying  acoustic  conditions  (such  as  thermal  structure, 
geographic  location,  time  of  day  and  year,  etc.)  it  should 
be  possible  to  determine  how  C(O)   and   At*   vary  as  a 
function  of  the  acoustic  environment. 

However,  it  is  not  necessary  to  observe  the  entire 
process  to  determine   C(0)   and   At  .   Specifically,  if 
the  process  of  transmission  loss  is  observed  at  two  points 
in  time  separated  by   A  t   and  these  pairs  of  points  are 
observed  many  times  at  varying  values  of   At   for  the  same 
process,  or  equivalently  from  processes  with  the  same 
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autoco variance  functions,  then  this  information  may  be  used 
to  estimate  C  (  A  t)   and  hence,   C(0)   and   At*.   This  may 
be  accomplished  in  the  following  way: 

Let 

Y(At)  =  X(t)  -  X(t+  At) 

represent  the  difference  between  two  measurements  in  trans- 
mission loss  taken  in  the  same  location  but  separated  in 
time  by   At. 

Then 

E(Y(  A  t))  =  0 

and 

Var(Y(  A  t) )  =  E(Y(  A  t)2) 

=  Var  (X(t))  +  Var  (X(t  +  At))  - 
2C  (At) 

=  2   (C(0)  -  C(  A  t)) 

The  value  of  Var(Y(  At))   may  then  be  determined  for  a 
given   At   from  the  observations  of  Y(  A  t)   corresponding 
to  pairs  of  observations  of  transmission  loss  with  the  same 
autocovariance  functions,  C(  A  t) .   C(0)   may  then  be 
represented  by 

C(0)  =  Var(Y(  A  t*))/2 
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where   Var  (Y(  At*))   is  the  value  of  Var  (Y( A  t) )   when 
C(  A  t)   approaches  zero  or  when  Var(Y(  At))   no  longer 
increases  with   At.   At*   is  then  the  time  interval  where 
this  first  occurs. 

If  pairs  of  observations  of  transmission  loss  cannot  be 
obtained  from  the  same  process,  then  it  is  necessary  to  pool 
information  from  different  processes  with  the  same  auto- 
covariance  function0   The  problem  here  arises  in  specifying 
the  range  of  conditions  under  which  this  may  be  done.   The 
approach  taken  in  this  study  was  to  characterize  the  local 
acoustic  environment  by  its  thermal  structure  where  it  was 
hypothesized  that  environments  with  similar  thermal  struc- 
tures will  produce  processes  with  similar  if  not  identical 
autocovariance  functions.   This  assumes  that  other  variables 
such  as  transmitter  and  receiver  depths  and  transmission 
ranges  will  remain  constant. 
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III.   NATURE  OF  THE  DATA 

The  data  utilized  in  this  report  was  made  available  by 
the  U.  S.  Underwater  Sound  Laboratory,  New  London, 
Connecticut,  and  reflects  the  raw  data  collected  during 
the  Acoustic,  Meteorological  and  Oceanographic  Survey  (AMOS) 
conducted  from  June,  1949,  through  April,  1953.   There  were 
nine  cruises  staged  during  these  four  years  which  covered 
the  North  Atlantic,  the  Norwegian  Sea,  and  to  a  lesser 
extent,  the  Mediterranean  Sea. 

Two  ships  were  employed  during  each  cruise  where  one 
acted  as  the  transmitting  platform  and  the  other  as  the 
receiving  platform.   Each  of  the  cruises  was  divided  into 
several  widely  spaced  stations  which  served  as  the  focal 
points  for  data  measurement  and  collection.   Within  each 
station,  the  acoustic  data  was  collected  at  several  trans- 
mission ranges  between  800  and  30,000  yards.   At  each 
transmission  range  the  transmission  loss  measurements  were 
observed  for  four  transmitting  frequencies  at  various 
source/receiver  depth  combinations.   The  transducer  depths 
varied  from  20  to  500  feet.   The  specific  data  recorded  was 
transmission  loss,  in  decibels  (db) ,  as  a  function  of 
source/receiver  depths  and  range  between  ships.   Table  X 
lists  the  information  that  was  recorded. 

Accuracy  of  the  data  becomes  particularly  important  in 
consideration  of  the  accuracy  of  any  results  obtained 
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through  the  use  of  that  data.   The  following  list  gives  an 
approximate  measure  of  the  accuracy  of  the  AMOS  data: 

1)  time  -  to  the  nearest  minute 

2)  range  -  to  the  nearest  5  yards 

3)  transducer  depths  -  to  the  nearest  foot 

4)  water  depth  -  to  the  nearest  10  fathoms 

5)  latitude  and  longitude  -  to  the  nearest  minute 

although  it  is  possible  that  the  two  ships 
may  have  drifted  as  much  as  a  mile  or  two 
from  the  recorded  position 

6)  BT  pattern  code  -  reflects  the  gross  thermal 

conditions  closest  to  the  time  that  the 
data  was  recorded. 
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IV.   ASSUMPTIONS 

As  indicated  in  the  introduction,  this  paper  deals  with 
determining  the  variation  of  transmission  loss  by  examining 
the  variation  of  the  change  in  transmission  loss  as  a 
function  of  the  elapsed  time  between  a  pair  of  observations. 
Ideally  it  would  be  best  if  this  change  in  decibel  loss 
could  be  measured  from  a  pair  of  observations  with  identi- 
cal source  and  identical  receiver  depths,  and  with  the  same 
horizontal  transmitting  range,  but  separated  by  an  interval 
of  time,  A  t.   Although  data  of  this  nature  was  not  available 
from  the  AMOS  experiment,  pairs  of  observations  of  db   loss 
were  obtained  having  the  following  characteristics: 

1)  One  of  the  observed  pair  measured  the 

transmission  loss  of  a  signal  generated 
at  depth  d-|_  and  received  at  depth  d2 
while  the  other  measured  transmission  loss 
of  a  signal  generated  at  depth  d2  and 
received  at  depth  d^. 

2)  Identical  transmitting  ranges  between  the 

source  and  receiving  ships  was  not  always 
available  within  a  pair  of  observations. 

Figures  (1)  and  (2)  depict  the  nature  of  an  ideal  pair 
of  observations  and  that  of  the  available  data. 

The  fact  that  the  change  in  transmission  loss  was 
measured  using  transmission  paths  that  'crossed'  instead  of 

15 


<>  Water 
Depth 

Source 


-1st  path 

"2nd  path 


Receiver 


Ship  1  ship  2' 

Figure  1.   Ideal  Transmission  Paths 


Ship 
Position 


Water 
Depth 


Source 


Receiver 


Alternative 
2nd  path 


Source 


Receiver 


Source 


Ship  1 


Ship  2    Ship  2   ship 
Position 


-  Figure  2.   Available  Transmission  Paths 

being  'parallel'  was  considered  not  to  be  a  serious  depart- 
ure from  the  ideal  conditions  (refs.  3  and  5).   This  follows 
from  the  assumption  of  a  generally  constant  horizontal 
thermal  environment,  and  therefore  an  acoustically 
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homogeneous  environment  over  the  range  of  operation  of  the 
two  ships  on  a  given  station,  and  theoretical  considera- 
tions which  postulate  that  transmission  loss  over  the  same 
path  from  different  directions  will  be  the  same.   It  is 
foreseen,  however,  that  this  approximation  to  the  ideal  case 
of  parallel  transmission  paths  will  increase  the  observed 
Var  (Y  (  A  t) )   for  small   At  due  to  a  contribution  caused 
by  small  scale  spacial  variability,  and  may  also  cause  a 
slight  underestimation  of   At.   However,  the  estimation 
of  C(0),  the  measure  of  total  variation  experienced  within 
a  generally  homogeneous  acoustic  environment,  will  not  be 
effected. 

The  second  data  characteristic  of  range  differences 
between  the  two  observations  in  a  pair  was  treated  in  the 
following  manner.   The  observed  change  in  transmission  loss 
was  adjusted  for  the  effect  of  this  range  differential  on 
the  basis  of  a  piecewise  linear  approximation  of  transmis- 
sion loss  vs.  horizontal  transmission  range.   This  linear 
approximation  seemed  reasonable  in  the  interest  of  keeping 
the  range  adjustment  computation  simple  and  yet  maintaining 
sufficient  accuracy.   Plotting  transmission  loss  vs.  trans- 
mission range  for  each  source/receiver  depth  combination  at 
each  station  revealed  that  the  resulting  relationships 
exhibited  the  same  general  shape  and  that  these  relation- 
ships could  be  approximated  by  two  linear  segments  with  the 
break  point  at  5,000  yards. 
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As  the  slopes  of  these  two  linear  segments  are  a  function 
of  the  local  environmental  conditions,  and  thus  varied  from 
station  to  station,  it  was  necessary  to  estimate  the  slope 
for  each  segment  at  each  station  for  each  source/receiver 
depth  combination.   This  was  done  on  the  basis  of  the  trans- 
mission loss  observed  at  the  various  ranges  occupied  on  a 
given  station.   Then  the  change  in  transmission  loss 
observed  at  a  particular  station  and  source/receiver  depth 
was  adjusted  using  the  corresponding  slope  of  transmission 
loss  vs.  transmission  range  in  the  following  manner: 

Adbadjusted  =  A  ^observed  "  b  X  A  ran^e 

where   A  db   represents  the  change  in  transmission  loss,   b 
is  the  slope  of  the  corresponding  linear  segment,  and 
A  range   is  the  range  differential  within  a  pair  of  obser- 
vations . 

In  some  cases  insufficient  data  was  available  to 
estimate  the  slope  of  one  or  both  segments  at  a  given 
station  and  source/receiver  depth  combination.   In  such 
cases  the  corresponding  adjustments  were  made  on  the  basis 
of  a  slope  which  was  typical  of  other  stations  on  the  same 
and  other  cruises.   Table  I  depicts  a  representative  cross- 
section  of  the  average  slopes  from  the  available  data. 
Listed  below  is  a  summary  of  the  typical  slope  values  selec- 
ted for  use  within  this  study  where  sufficient  data  was  not 
available  for  slope  calculations. 
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Applicable 
Range  Interval      Slope  Value      Cruises 

Above  .0015  5,  7,  8,  9,  10, 

5,000  yards  11,  12 

Below  .0045  7,  8,  9,  10,  12 

5,000  yards         .0100  11 

.0280  5 


Admittedly  this  adjustment  procedure  is  subject  to 
error.   In  particular,  the  error  incurred  by  using  a  linear 
approximation  in  the  first  place,  and  secondly  in  the  esti- 
mation of  the  slope  of  the  transmission  loss  vs.  transmission 
range  relationship.   It  was  felt,  however,  that  in  general 
the  total  error  inherent  in  this  procedure  could  be 
restricted  to  at  most  ll.5db  by  making  no  adjustments  with 
range  differentials  within  pairs  of  observations  of  more 
than  500  yards.   Errors  of  such  magnitude  were  considered 
acceptable  in  light  of  the  original  accuracy  in  measuring 
transmission  loss  in  the  AMOS  experiment.   In  actual 
practice,  this  error  was  further  limited  by  the  fact  that 
few  adjustments  in  change  in  transmission  loss  were  made  for 
range  differentials  greater  than  250  yards. 
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V.   RESULTS 

A.   METHODOLOGY  OF  INVESTIGATION 

The  purpose  of  this  research  is  to  investigate  the 
nature  and  magnitude  of  the  variability  in  sound  propogation 
loss  as  characterized  by   C (  A  t)   by  investigating  the 
variability  in  the  change  in  transmission  loss  as  a  function 
of  time  between  observations.   The  calculation  of  trans- 
mission loss  as  a  function  of  time,  denoted  by  Y(  A  t), 
was  obtained  from  the  AMOS  data  in  all  instances  where  a 
pair  of  transmission  loss  observations  were  taken  under  the 
following  conditions: 

1)  The  compared  observations  were  taken  at 

the  same  station. 

2)  The  range  differential  between  observations 

was  no  greater  than  500  yards. 

3)  The  source  depth  of  one  transmission  equalled 

the  receiver  depth  of  the  other,  and 
vice-versa,  thus  resulting  in  a  compari- 
son of  observations  with  similar  but  not 
identical  transmission  paths. 

These  calculations  of  change  in  transmission  loss  over 
common  paths  were  indexed  by  the  following  environmental 
and  non-environmental  factors: 

1)  Source  and  receiver  depths 

2)  Transmission  range 
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3)  Transmitting  frequency 

4)  Water  depth 

5)  Mixed  layer  depth 

6)  Geographical  location 

7)  Time  interval  between  corresponding  observations 

8)  Time  of  day  and  year 

The  investigation  of  variability  in  sound  transmission 
was  confined  to  the  influence  of  the  above  factors.   To 
determine  this  influence  using  the  available  data,  it  was 
necessary  to  pool  data  collected  under  varying  environ- 
mental conditions.   Such  pooling  is  only  applicable,  however, 
over  a  range  of  environmental  conditions  where  the  acoustic 
variability  is  characterized  by  a  common  autocovariance 
function.   Hence,  it  is  necessary  to  determine  how  to  iden- 
tify the  range  of  conditions  which  lead  to  a  common  auto- 
covariance function.   In  particular,  it  is  necessary  to 
determine  which  factors,  both  environmental  and  non- 
environmental,  influence  the  nature  and  magnitude  of  the 
variation  of  propogation  loss  in  the  ocean  and  then  pool 
over  the  entire  range  of  those  factors  which  did  not  so 
influence  this  variation. 

The  procedure  followed  was  to  plot   Var(Y(  At))   vs. 
At   using  the  entire  sample  set  of  data  and  thus  represent- 
ing the  entire  range  of  environmental  and  non-environmental 
conditions  encountered.   This  entire  sample  of  data  was 
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then  subdivided  into  sub-samples  on  the  basis  of  the 
source/receiver  depth  combination  and   Var(Y(  A  t) )   vs. 
At  was  plotted  for  each  sub-sample.   By  this  procedure, 
it  was  determined  that  the  nature  and  magnitude  of 
Var(Y(  At))   depended  upon  the  sub-sample  (i.e.,  source/ 
receiver  depth  combination)  considered.   Hence,  it  was 
deemed  inappropriate  to  pool  the  data  over  source/receiver 
depth  combinations  when  considering  the  influence  of  other 
factors  on   Var(Y(A  t) )  . 

To  determine  if  the  thermal  structure  influenced 
Var(Y(  At)),  each  of  the  above  sub-samples  was  further 
subdivided  on  the  basis  of  the  location  of  the  receiver  and 
source  depths  with  respect  to  mixed  layer  depth.   On  the 
basis  of  this  subdivision,  it  was  determined  that  the  posi- 
tioning of  the  source  and  the  receiver  with  respect  to  the 
mixed  layer  depth  also  influenced  the  nature  and  magnitude 
of   Var(Y(  At)). 

The  influence  of  geography  and  transmission  range  were 
also  investigated  in  the  same  manner  by  subdividing  the 
sample  data  on  the  basis  of  these  factors.   As  it  was 
necessary  to  keep  the  respective  sub-samples  as  large  as 
possible  when  considering  a  given  factor,  it  was  found  that 
certain  sub-samples  should  not  be  further  subdivided.   Such 
was  the  case  for  some  of  the  less  representative  source/ 
receiver  depths. 


22 


By  proceeding  in  this  manner,  it  was  determined  that,  of 
the  factors  being  considered,  those  which  were  judged  to 
have  a  significant  influence  on  the  nature  and  magnitude  of 
Var(Y(  At))   were:   1)  source/receiver  depth,   2)  location 
of  the  source  and  receiver  with  respect  to  the  mixed  layer, 
3)  geographic  location,  and  4)  transmission  range.   It  was 
also  determined  that  the  transmitting  frequency  had  a  marked 
effect  on  the  magnitude  of  Var(Y(  At)).   The  actual  nature 
of  these  influences  will  be  discussed  later. 

It  should  be  emphasized  that   Var(Y(  At))   is  twice 
the  value  of  C(0)   and  thus   Var(Y(  t) )  ,   the  variation  of 
transmission  loss. 

B.   LIMITATIONS  EFFECTING  RESULTS 

There  are  two  considerations  which  tended  to  limit  the 
results  of  this  paper.   The  first  was  that  the  amount  of 
data  available  to  estimate   Var(Y(  At))   as  a  function  of 
At  was  limited  in  some  of  the  finer  subsamples  defined  in 
the  previous  section.   Hence  it  was  found  necessary,  in 
order  to  obtain  a  reasonably  precise  estimate  of 
Var(Y(  At)),  to  pool  the  available  data  into  several  con- 
secutive 15  minute  intervals. 

The  second  consideration  was  that,  due  to  physical 
limitations  imposed  on  observation  of  the  data  during  the 
AMOS  experiment  for  a  given  source/receiver  depth  combina- 
tion, the  majority  of  the  available  data  was  restricted  to 
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at  most  three  15  minute  intervals.   These  three  intervals, 
however,  were  not  necessarily  the  same  for  each  source/ 
receiver  depth  combination,  but  tended  to  the  small  values 
of  At   for  closely  spaced  source  and  receiver  depths  and 
to  the  large  values  for  widely  spaced  source  and  receiver 
depths. 

The  implication  of  these  data  limitations  was  that  A  t 
could  not  be  determined  with  a  precision  greater  than  15 
minutes o   However,  this  was  thought  not  to  be  a  serious 
limitation.   A  possibly  more  serious  limitation  is  that, 
because  of  the  dispersion  of  the  data  as  a  function  of   At, 
the  value  of   Var(Y(  A  t) )   cannot  be  estimated  with  equal 
precision  from  interval-to-interval  and  from  depth 
combination-to-depth  combination.   In  particular,  for  a 
given  source/receiver  depth  combination,  there  may  be  little 
or  no  data  available  for  the  interval  when  Var(Y(  A  t) ) 
reaches  its  maximum  value.   In  general  it  is  felt  that  this 
did  not  happen,  although  there  is  no  way  of  verifying  this 
at  the  present  time. 

C.   FINDINGS 

1.   Distribution  of  Transmission  Loss 

As  suggested  in  the  introduction,  knowledge  of  the 
nature  of  the  variability  in  transmission  loss  may  be 
utilized  to  produce  a  realistic  simulation  of  the  acoustic 
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environment.   To  do  this  it  is  necessary  to  know  not  only 
the  magnitude  of  this  variation,  but  also  its  distribution. 

The  typical  assumption  that  this  distribution  is  normal 
is  evaluated  on  the  basis  of  the  available  data  for  a 
representative  cross-section  of  subsamples  of  the  data. 
The  basis  for  this  examination  is  the  Kolomogorov-Smirnov 
one-sample  goodness-of-f it  test  (ref .  2) .   Table  II  gives 
the  results  for  all  subsamples  considered,  and  the  range  of 
conditions  which  define  the  tested  subsample.   Also  given 
are   Dmax,   the  maximum  absolute  difference  between  sample 
cumulative  distribution  and  a  normal  cumulative  distribution 
with  mean  zero  and  sample  variance  of  the  subsample.   The 
critical  value  has  the  same  connotation  usually  associated 
with  hypothesis  testing  and  represents  the  level  of  signifi- 
cance at  which  the  null  hypothesis  of  no  departure  from 
normality  would  be  rejected. 

The  general  conclusion  on  the  basis  of  these  results  is 
that  the  normal  distribution  is  an  adequate  approximation 
for  the  distribution  of  transmission  loss. 

2.   Factors  Effecting  the  Variation  of  Transmission  Loss 

a.   Source/Receiver  Depth 

The  first  attempt  at  determining  which  environ- 
mental and  non-environmental  factors  influenced  the  nature 
and  magnitude  of   Var(Y(  A  t) )   was  in  the  area  of 
source/receiver  depths.   The  available  data  was  subdivided 
on  the  basis  of  several  source/receiver  depth  combinations. 
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The  results  thus  obtained  proved  to  be  quite  noisy  and  it 
was  decided  to  employ  the  same  procedure  for  cruise  10  only 
since  the  thermal  structure  appeared  more  consistent  on  that 
cruise.   The  results  of  this  procedure  seemed  satisfactory 
for  exposition  purposes  and  are  summarized  in  Table  III. 

It  is  interesting  to  note  that  not  only  does   Var(Y(At)) 
increase  as  time  increases,  but  that  this  variation  also 
increases  with  an  increasing  depth  differential  between 
source  and  receiver.   Also  noted  was  that   Var(Y(  A  t) ) 
tended  to  decrease  as  the  source  depth  increased.   The 
latter  is  probably  attributable  to  surface  effects, 
b.   Mixed  Layer  Depth 

Once  it  had  been  decided  that  the  different 
source/receiver  depth  combinations  each  influenced 
Var(Y(  A  t) )   differently,  the  next  step  was  to  consider 
the  position  of  the  source  and  receiver  with  respect  to  the 
depth  of  the  mixed  layer.   This  was  accomplished  by  further 
subdividing  the  data  into  two  subsamples  on  the  basis  of 
whether  both  the  source  and  receiver  were  in  the  mixed 
layer,  or  at  least  one  was  below  the  mixed  layer.   The 
rational  behind  this  subdivision  was  to  determine  if  the 
variability  of  transmission  loss  was  different  when  the 
transmission  path  was  above  the  thermocline  or  across  the 
thermocline. 

Table  IV  shows  the  results  of  this  investiga- 
tion and  clearly  indicates  that   Var(Y(At))   is 
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considerably  less  within  the  mixed  layer.   The  differences 
observed  may  in  part  be  a  function  the  existence  of  internal 
waves. 

c.   Geographic  Location 

Geographic  location  was  found  to  have  a 
significant  effect  on  the  variation  of  transmission  loss. 
To  facilitate  investigating  the  nature  of  this  effect,  areas 
were  selected  which  were  thought  to  have  approximately  the 
same  large  scale  environmental  characteristics.   The  areas 
chosen  were: 

1)  50-80  N,  all  longitudes 

2)  0-50  N,   0-60  W 

3)  0-50  N,  60-90  W 

Four  source/receiver  depth  combinations  were  then 
picked  which  would  ensure  that  sufficient  pairs   of  observa- 
tions would  be  available  within  each  of  these  depth  combina- 
tions.  The  subsamples  have  been  further  subdivided  accord-t- 
ing to  the  position  of  the  source  and  receiver  with  respect 
to  the  mixed  layer.   Table  V  depicts   Var(Y(  A  t) )   for  the 
geographic  areas  and  the  four  depth  combinations  selected. 
The  variation  of  transmission  loss  for  each  combination 
pooled  over  all  latitudes  and  longitudes  has  also  been 
included  in  the  table  for  comparative  purposes. 

As  indicated  by  the  entries  in  Table  V,  there 
does  not  appear  to  be  any  consistently  predictable  direc- 
tion or  magnitude  associated  with   Var(Y(  At))  across  the 
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three  geographic  locations.   However,  the  variations  do  seem 
to  change  significantly  both  from  the  "all  latitude- 
longitude"  case,  and  across  geographic  areas.   This  would 
seem  to  substantiate  the  hypothesis  that  geographic  loca- 
tion is  a  significant  factor  effecting   Var(Y(  At)).   The 
data  in  Table  VI,  which  breaks  the  variation  of  trans- 
mission loss  down  into  15  minute  intervals,  also  appears 
to  substantiate  this  contention. 

d.  Transmitting  Frequency 

The  effect  of  transmitting  frequency  on 
Var(Y(  At))   is  apparent;  increasing  the  frequency  in- 
creased the  magnitude  of  the  variation.   Two  frequencies 
were  selected  to  investigate  this  hypothesis,  8  and  25  KC. 
Table  VII  is  a  summary  of   Var(Y(  At))   for  both  frequen- 
cies for  several  source/receiver  depth-geographic  location- 
mixed  layer  depth  combinations.   The  increase  in  Var (Y ( A  t) ) 
is  evident. 

e.  Transmission  Range 

The  transmission  range  appeared  to  have  a 
significant  effect  on   Var(Y(  At)).   Table  VIII  shows  this 
variation  for  range  intervals  of  0-10  Kyds  and  10-30  Kyds 
for  pairs  of  observations  at  the   8  KC  transmitting  frequency 
pooled  appropriately  according  to  source/receiver  depth, 
geographic  location  and  mixed  layer  depth.   An  examination 
of  this  data  suggests  that  Var(Y(  A  t) )   was  generally 
lower  for  greater  ranges  and  that  estimates  of  variation 
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in  the  mixed  layer  was  less  consistent  than  for  those 
transmissions  crossing  the  layer. 

3.   Investigation  of  Large  Values  of  Change  in 
Transmission  Loss 


Because  of  the  limited  amount  of  data  in  some 
subsamples  considered,  the  estimator  of  Var(Y(  A  t) )   was 
quite  sensitive  to  a  not  uncommon  occurrence  of  dispropor- 
tionately large  values  of  the  change  in  transmission  loss. 
In  an  attempt  to  determine  the  nature  of  the  circumstances 
which  lead  to  such  observations,  a  sample  of  observed 
changes  in  transmission  loss  whose  absolute  value  exceeded 
9db  were  examined  as  to  the  conditions  under  which  they  were 
observed.   This  sample  consisted  of  all  instances  for 
20/100,  50/100,  50/250  and  100/500  foot  source/receiver 
depths  at  8  KC  transmitting  frequency,  and  20/100  and  100/500 
foot  source/receiver  depths  at  25  KC  transmitting  frequency. 
Table  IX  is  a  summary  of  these  data. 

On  examination  of  these  data  it  is  concluded  that, 
except  for  a  marked  effect  of  frequency  and  a  suggested 
effect  of  time  of  day,  no  common  single  or  set  of  circum- 
stances is  related  to  the  large  changes  in  transmission 
loss.   This  does  not  mean  that  such  a  set  of  circumstances 
does  not  exist,  but  that  it  is  not  apparent  what  such  a  set 
would  be  based  on  in  consideration  of  these  data  and  this 
specific  set  of  indices. 
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This  latter  point  has  lead  to  the  consideration 
that  there  may  be  a  more  appropriate  set  of  indices  with 
which  a  measured  change  in  transmission  loss  may  be  classi- 
fied.  The  nature  of  such  a  set  of  indices  is  suggested  by 
considering  the  results  of  the  present  study.   That  is,  that 
Var(Y(  At))   varied  with  the  proximity  of  the  source  or  the 
receiver  to  the  surface,  the  depth  differential  between  the 
source  and  receiver,  the  relationship  of  the  source  and 
receiver  to  the  mixed  layer,  geographic  location,  etc. 
These  factors  are  exactly  those  which  influence  the  mode  of 
transmission  of  a  signal  in  the  general  transmission  loss 
equations  obtained  in  the  AMOS  study  (ref.  5). 

In  that  study  it  was  determined  that  transmission 
loss  depended  on  the  mode  of  transmission.  In  particular 
on  whether  transmission  was  by 

1)  Direct  path 

2)  Single  surface  bounce 

3)  Multiple  surface  bounce 

4)  Leakage  by  diffraction  or  scattering 

5)  Depressed  sound  channel 

6)  Bottom  bounce 

On  the  basis  of  the  results  of  the  present  study, 
it  appears  that  the  nature  and  magnitude  of  the  variability 
in  transmission  loss  might  also  vary  with  transmission  mode. 
Further,  it  is  envisioned  that  subdividing  the  available 
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data  on  the  basis  of  transmission  mode  will  yield  a  more 
uniform  distribution  of  At   for  a  given  subsample  and,  in 
general,  increase  the  amount  of  data  available  for  the  esti- 
mation of  Var(y(  At))   for  any  particular  interval  of  At. 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  aim  of  this  study  was  to  determine  the  environmental 
and  non-environmental  factors  which  influence  the  nature 
and  magnitude  of  the  variability  in  transmission  loss,  and 
then  to  characterize  this  variability  as  a  function  of  these 
influencing  factors. 

On  the  basis  of  this  study  it  was  found  that  the 
variation  in  transmission  loss  was  influenced  by  the  depth 
of  the  source  and  receiver  of  the  transmitted  signal,  the 
relationship  of  the  source  and  receiver  depths  with  the 
mixed  layer  depth,  geographical  location,  and  transmission 
range. 

The  nature  of  this  influence  was  that  the  magnitude  of 
variability  increased  with  increasing  transmitting  frequency, 
with  increasing  depth  differential  between  the  source  and 
receiver,  with  the  increasing  closeness  of  the  source  and 
receiver  to  the  air-sea  interface,  with  decreasing  trans- 
mission range,  and  when  transmitting  across  the  thermocline 
as  opposed  to  transmitting  exclusively  within  the  mixed 
layer. 

The  time  interval  between  uncorrelated  observations  was 
found  to  be  shortest  near  the  surface  as  opposed  to  at 
depth,  to  increase  with  increasing  source/receiver  depth 
differential,  and  to  decrease  when  transmitting  across  the 
thermocline  as  opposed  to  transmitting  entirely  within  the 

mixed  layer. 
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The  magnitude  of  the  variability  in  transmission  loss 

experienced  under  the  varying  conditions  ranges  over  one 

2 
order  of  magnitude  from  10  db   in  the  near  horizontal  trans- 

2 

mission  in  the  mixed  layer  at  8  KC  to  100  db   crossing  the 

thermocline  in  the  25  KC  case.   Such  figures  are  only  esti- 
mates and  as  such  are  subject  to  statistical  variation,  but 
the  trends  in  magnitude  of  the  variability  appear  to  be 
consistent  in  the  various  cases  considered.   Therefore,  the 
magnitudes  of  the  variability  determined  in  this  study  may 
be  used  as  approximations  of  the  true  variability  which  will 
be  experienced  in  the  corresponding  environmental  situation. 

In  conclusion,  it  is  observed  that  the  factors  which 
influenced  the  nature  and  magnitude  of  the  variability  in 
transmission  loss  are  exactly  the  same  as  those  found  in  the 
AMOS  study  to  influence  the  mean  transmission  loss  (ref .  5) . 
In  particular,  the  factors  of  source  and  receiver  depth, 
mixed  layer  depth,  transmission  range  and  surface  tempera- 
ture were  used  to  determine  the  mode  of  transmission  of  a 
given  transmitted  signal.   The  transmission  loss  expected 
was  then  calculated  using  an  equation  peculiar  to  this  speci- 
fic transmission  mode. 

In  a  li>:e  manner  it  was  suggested  that  the  nature  of  the 
variability  in  transmission  loss  may  be  more  adequately 
classified  and  the  magnitude  of  the  variability  more  pre- 
cisely specified  if  the  observed  change  in  transmission 
losses  were  further  indexed  and  examined  by  transmission 

33 


mode.   It  is  therefore  recommended  that  this  be  the  next 
area  of  investigation  in  the  study  of  the  nature  and  magni- 
tude of  the  variability  of  transmission  loss  in  the  ocean. 
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APPENDIX  A 

I.   DESCRIPTION  OF  DATA 

A.   TAPE  RECONSTRUCTION 

The  information  as  received  from  the  U.  S.  Underwater 
Sound  Laboratory  was  in  the  form  of  9,  seven-track  tapes 
containing  appropriate  spacial,  temporal,  and  acoustic 
information.   These  tapes  have  now  been  consolidated  into 
two  nine-track  tapes  suitable  for  mounting  on  the  standard 
IBM  tape  drive  units.   One  tape  contains  only  the  bathy- 
thermograph (BT)  information  while  the  other  contains  all 
other  necessary  data  including  acoustic  measurements 
pertinent  to  this  report.   Table  X  displays  the  information 
that  is  recorded  on  this  latter  tape.   The  following  informa- 
tion has  been  omitted  from  each  record  on  the  revised 
acoustic  tape,  but  is  still  available  on  the  original  seven- 
track  tapes: 

1)  Ambient  noise  indicator 

2)  Wind  force  (beaufort) 

3)  Unknown  contents  following  water  depth  entry 
The  acoustic  tape  contains  data  for  cruises  5,  7,  8,  9, 

10,  11,  and  12.   Each  record  is  7  2  characters  in  length 
although  each  only  utilizes  60  character  spaces.   The  addi- 
tional 12  spaces  have  been  left  blank  to  allow  room  for  any 
additional  information  deemed  appropriate  and  to  facilitate 
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placing  this  data  on  a  user  disc  if  desired.   No  special 
records  have  been  added  to  indicate  the  end  of  a  cruise  or 
station.   It  should  be  noted  that  all  field  entries  are  in 
terms  of  integer  numbers  with  the  exception  of  fields  5, 
24,  and  27  which  will  require  character  format  (A-format 
under  FORTRAN  IV) . 

The  data  for  cruises  2  and  3  were  deemed  inappropriate 
and  therefore  not  included  on  the  acoustic  tape.   The 
reason  was  that  the  source/receiver  depths  and  two  of  the 
four  transmitting  frequencies  do  not  correspond  to  the 
requirements  of  this  study.   That  is,  both  cruises  show 
source  depths  set  at  a  constant  15  feet  while  the  receiver 
depth  is  varied.   Consequently,  none  of  this  data  would  be 
accepted  by  the  computer  program.   The  program  requires  the 
comparison  of  one  record  with  source  depth  at  d.  feet  and 
the  receiver  depth  of  d2  feet,  with  another  record  having 
the  source  at  depth  d2  and  receiver  at  depth  d-j_. 

B.   COMBINING  ACOUSTIC  DATA 

An  examination  of  the  original  tapes  revealed  that 
either  1,  2,  or  3  acoustic  measurements  were  taken  for  each 
particular  location,  time,  and  range.   Hence,  some  method 
must  be  employed  to  combine  these  sets  of  data  in  order  that 
they  be  consistent.   One  approach  would  have  been  to  take  an 
average  each  time  a  grouping  of  two  or  three  measurements 
was  encountered.   Another  approach,  and  the  one  subsequently 
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adopted,  would  be  to  simply  sum  the  acoustic  data  for 
groupings  of  three  measurements,  sum  and  multiply  by  3/2 
for  groupings  of  two,  and  triple  any  single  entries.   The 
resulting  figures  would  then  be  entered  on  the  revised 
acoustic  tape  in  place  of  the  old  transmitting  frequency 
data,  giving  a  single  entry  for  each  particular  location, 
time,  and  range.   This  latter,  method  was  utilized  in  the 
interest  of  conserving  valuable  computer  time,  and  in 
consideration  of  computational  and  programming  ease.   In 
summary,  all  acoustic  field  entries  at  each  of  the  four 
transmitting  frequencies  must  be  divided  by  3  before  they 
can  be  used  in  any  computations. 

C.   BATHYTHERMOGRAPH  PATTERN  CODE  AND  ITS  USE 

The  bathythermograph  (BT)  pattern  code  is  a  two  digit 
number  and  will  be  found  in  fields  18  and  20  on  the  acoustic 
tape. 

Interpretation  of  the  code  has  been  simplified  for  use 
within  this  study.   Figure  (3)  indicates  the  general  nature 
of  the  thermal  structure.   To  use  the  figure,  enter  with 
the  first  digit  of  the  code  across  the  top.   This  gives  a 
picture  of  the  approximate  temperature  gradient  down  to  100 
feet  (associated  depths  may  be  read  in  the  '1st  digit* 
column).   Next,  enter  with  the  second  digit.   This  gives  a 
picture  of  the  approximate  temperature  gradient  from  100  to 
250  feet.   Connecting  the  two  pictures  thus  obtained  results 
in  a  complete  recap  of  the  thermal  conditions  from  0  to  250 
feet. 
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Figure  3.  -  BT  Pattern  Code  Interpretation 

This  BT  pattern  code  information  was  used  to  establish 
the  position  of  the  mixed  layer. 
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APPENDIX  B 

I.   EXPLANATION  OF  COMPUTER  PROGRAM 
This  section  has  been  included  in  the  report  in  order 
that  the  program  may  be  used  as  a  base  for  further  studies 
in  the  field  of  temporal  and  spacial  effects  on  the  varia- 
tion of  decibel  loss  underwater.   The  entire  program  has 
been  written  and  compiled  under  FORTRAN  IV,  G-level. 

A.   GENERAL  EXPLANATION 

The  program  is  designed  to  read  one  data  record  at  a 
time  from  the  acoustic  tape,  check  the  appropriate  data 
against  the  parameters  established  by  the  input  array,  and, 
if  possible,  pass  the  entire  record  to  subroutine  C0MP1. 
Subroutine  C0MP1  then  separates  the  records  passed  to  it 
into  four  categories  according  to  range  between  the  trans- 
mitting and  receiving  ships,  the  range  break  point,  and 
source  and  receiver  depths.   For  example,  if  the  source/ 
receiver  depths  are  either  20/100  or  100/20  feet  and  the 
range  break  point  is  5,000  yards,  the  categories  are: 

Range  Interval 

-   5,000  yards 

^-5,000  yards 

5  5,000  yards 

y  5,000  yards 
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Source/Receiver 

Category 

Depth 

I 

20/100 

II 

20/100 

III 

100/20 

IV 

100/20 

Once  this  data  has  been  collected  for  a  single  station,  the 
appropriate  slope  information  is  computed  by  subroutine 
RGRESS  and  the  next  station  is  then  processed  through  the 
main  program,  C0MP1,  and  RGRESS. 

At  the  end  of  each  cruise,  the  processed  records  are 
sent  to  subroutine  C0MPo   This  routine  computes  the  change 
in  transmission  (db)  loss  in  a  sequential  manner.   That  is, 
a  record  can  never  be  used  in  a  computation  more  than  once. 
The  computation  of   A  db   loss  proceeds  in  the  following 
manner: 

1)  At  any  given  time,  two  adjacent  records  are 

considered. 

2)  The  program  computes  the  difference  between 

ranges  of  the  two  records.   If  this  difference 
is  greater  than  the  specified  amount  '  Z1,  then 
no  computation  is  made  and  the  next  data  set 
is  picked  up, 

3)  If  the  range  difference  is  less  than  or  equal 

to  the  specified  amount,  Z,  then  the  change 
in  db  loss  is  calculated. 

Of  the. two  records  being  considered  at  any  particular  time, 

the  decibel  loss  of  the  second  is  adjusted  in  the  direction 

of  the  first  record  with  respect  to  the  difference  in  range 

associated  with  the  two  records „   This  is  the  point  at  which 

the  slope  information  passed  from  subroutine  COMPl  is  util- 
ized.  The  only  exception  to  this  computational  procedure 
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is  when  one  range  falls  below  the  break  point  and  the  other 
above.   In  that  event,  both  db  losses  are  adjusted  to  the 
range  break  point.   Once  a  computation  has  been  completed, 
the  change  in  decibel  loss  and  the  time  difference  between 
the  two  records  are  stored  in  appropriate  arrays.   At  the 
end  of  each  cruise,  this  information  may  be  plotted  at  the 
users  option.   Also  within  this  subroutine  are  two  arrays 
which  collect  the  same  information  for  all  cruises  to 
facilitate  a  cumulative  plot. 

At  the  end  of  the  tape,  the  program  calls  subroutines 
SIGMA  and  BTLAYR.   SIGMA  simply  sorts  the   A db   losses  with 
respect  to  their  associated  time  increments  into  5  and  then 
15  minute  intervals,  computes  the  variance  within  each 
interval,  and  plots  the  variance  vs.  the  change  in  time. 
It  also  calculates  the  variance  and  standard  deviation  over 
all  data  passed  to  it  at  the  end  of  each  tape  iteration. 

BTLAYR  is  designed  to  consider  all  comparisons  which 
lead  to  computation  of   Adb   loss  in  subroutine  COMP,  and 
separate  the  pairs  of  observations  into  one  of  three  cate- 
gories:  1)  above  the  layer,   2)  transmitter  and  receiver 
both  below  the  layer,  and   3)  all  transmissions  with  the 
receiver  and  source  neither  above  or  below  the  thermocline. 
Once  all  records  have  been  categorized,  the  change  in  db 
loss  vs.  change  in  time  may  be  plotted  at  the  users  option. 
Next  subroutine  SIGMA  is  called,  which  computes  the 
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appropriate  output  as  described  above  for  each  of  the  three 
transmission  categories. 

Subroutine  SIME  solves  simultaneous  linear  equations 
necessary  to  subroutine  RGRESS.   OSPLOT  is  a  standard  IBM 
plotting  package  and  is  used  exclusively  for  plotting  within 
the  program. 

B.  INPUT  SPECIFICATIONS  AND  DESCRIPTION 

An  input  array  has  been  established  in  order  to  lend  a 
certain  amount  of  flexibility  in  filtering  data  through 
specific  parameters.   This  information  must  be  punched  on 
computer  cards  for  inclusion  at  the  end  of  the  program  and 
will  be  either  integer  or  alphabetic  in  nature.   There  must 
be  six  entries  on  each  card  with  the  last  letter  of  the 
numeral  entered  successively  in  columns  6,  12,  18,  24,  30, 
36,  42,  48,  54,  60,  and  66.   The  last  9  entries  on  the  last 
card  may  be  left  blank.   Table  XI  fully  explains  what 
information  is  required  and  gives  a  general  description  of 
the  input  variables.   The  input  data  is  read  into  the  NK 
array.   Note  that  NK(15)  and  NK(17)  require  a  character 
format.   Tables  XII  and  XIII  describe  the  arrays  and  vari- 
ables \ised  within  the  computer  program. 

C.  READING  DATA  FROM  THE  ACOUSTIC  TAPE 

The  acoustic  tape  data  format  is  described  by  Table  X. 
There  is  one  peculiarity  which  exists  within  the  formatting 
structure.   Note  that  the  range  between  the  two  ships  is 
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divided  into  two  separate  field  entries.  In  this  case  it 
is  necessary  to  compute  the  range  in  the  following  manner 
(tape  records  are  loaded  into  the  IN  array) : 

Range  =  1N(9)  +  1,000  X  IN(8) 

All  data  on  the  tape  is  integer  in  nature  except  fields  5, 
6,  24  and  27,  which  are  alphameric. 

D.   ESTABLISHING  THE  MIXED  LAYER  DEPTH 

Within  each  pair  of  compared  data  sets  are  four  BT 
pattern  codes,  one  per  ship  per  record.   Since  there  exists 
a  change  in  time  and  a  range  differential  between  data  sets, 
it  seemed  appropriate  to  combine  the  depths  associated  with 
the  pattern  codes  in  some  manner.   One  acceptable  method  and 
the  one  adopted,  takes  an  average  of  the  associated  thermal 
layer  depths  across  records  first  for  the  transmitting  ship, 
and  then  the  receiving  ship.   This  averaging  is  done  in  the 
same  manner  for  determining  both  the  top  and  the  bottom  of 
the  thermocline.   Once  this  has  been  accomplished,  the  mini- 
mum of  the  two  top  of  the  layer,  and  the  maximum  of  the  two 
bottom  of  the  layer  averages  are  selected  to  represent  the 
upper  and  lower  bounds  of  the  thermocline  respectively. 

As  an  example,  consider  some  BT  pattern  code   X]^' 
where  it  was  desired  to  determine  both  the  top  and  the 
bottom  of  the  thermocline.   Computing  the  top  and  bottom  of 
the  thermocline  was  done  in  the  following  manner: 
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1)  Determining  the  top  of  the  thermocline  (enter 
with  X]_) 


h. 

Depth  (ft) 

*2 

Depth  (ft) 

i 

0 

1 

100 

2 

0 

2 

100 

3 

see  X2 

3 

100 

4 

0 

4 

100 

5 

30 

5 

150 

6 

30 

6 

150 

7 

50 

7 

200 

8 

see  X9 

8 

600 

2)  Determining  the  bottom  of  the  thermoline 
(enter  with  X2) 


ik 

Depth  (ft) 

S, 

Depth  (ft) 

i 

300 

1 

100 

2 

200 

2 

50 

3 

150 

3 

30 

4 

300 

4 

100 

5 

300 

5 

50 

6 

300 

6 

100 

7 

300 

7 

100 

8 

see  Xt 

8 

600 

Note  that  when  checking  for  the  top  of  the  thermocline  that 
when   Xj_  =  3  or  8,  this  implies  looking  at  the  second  digit. 
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Likewise,  when  determining  the  bottom  of  the  thermocline 
and  X2  =  8,  then  look  at  the  first  digit. 
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TABLE  I. 


Average  Slopes  Pooled  Over  All  Possible  Data 


Slope 


Source/Receiver  Cruise 
Depth  (ft.)      No. 


50/200 


30/150 


20/100 


50/400 


100/400 


ch/d. 


d9/d- 


-  5  Kyds  ^5  Kyds   -  5  Kyds  ^5  Kyds 


5 

.003113 

.06622 

.001893 

7 

.001978 

8 

.003122 

.001446 

9 

.00447  6 

.001040 

10 

.003580 

.001894 

11 

.010430 

.000683 

9 

.004950 

.000759 

.000358 

10 

.004975 

.001400 

.004202 

.001356 

5 

.028020 

.002022 

.020260 

.003018 

7 

O0019581 

.002236 

8 

.005097 

.001348 

O004221 

O015560 

9 

.009709 

.001062 

.008700 

.001550 

11 

.009000 

.001000 

.010100 

.001000 

8 

.003831 

.001400 

.003430 

.001358 

9 

.004455 

.001191 

10 

.003607 

.001924 

11 

.011110 

.000997 

5 

.000599 

7 

.000930 

8 

.004718 

.001453 

.001964 

.001313 

9 

.006369 

.001229 

10 

.003980 

.001434 

11 

.016850 

.001105 

:Source/receiver  depths 
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TABLE  IV. 


The  Effect  of  Mixed  Layer  on  Var  (y  (  /\   t) ) 


Source/Receiver 
Depth  (ft.) 

Above 
14.8 

** 

Crossing 

20/100 

28.4 

20/500 

57.2 

58.7 

30/150 

19.2 

- 

50/100 

20.5 

24.9 

50/250 

5.6 

25.7 

50/400 

18.6 

28.5 

50/500 

14.6 

34.3 

100/400 

11.7 

25.4 

100/500 

12.0 

41.5 

Both  source  and  receiver  in  the  mixed  layer 

** 

At  least  source  or  receiver  below  the  mixed  layer,  i.e., 

crossing  the  thermocline 
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TABLE  V. 

The  Gross  Effect  of  Geographic  Location 
on  Var  (  A  db  loss)  -  8  Kc 


Source/Receiver   Latitude   Longitude 
Depth  (ft.) 

20/100 


50/100 


50/250 


100/500 


(degs) 

(degs) 

*7C 

Crossing 

Above*** 

All 
50-80N 
0-50N 
0-50N 

All 
All 
0-60W 
60-90W 

28.30(192)* 
23.12(111) 
37.60(65) 
8.30(5) 

14.82(101) 

2.60(2) 
12.29(57) 
19.16(22) 

All 
50-80N 
0-50N 
0-50N 

All 
All 
0-60W 
60-90W 

24.95(216) 
20.00(115) 
29.96(79) 
25.95(11) 

20.53(167) 

6.77  (55) 
27.56(63) 
29.10(29) 

All 
50-80N 
0-50N 
0-50N 

All 
All 
0-60W 
60-90W 

25.71(151) 
23.54(94) 
20.52(35) 
43.26(22) 

7.53(48) 
4.86(39) 
24.61(6) 
7.89(3) 

All 
50-80N 
0-50N 
0-50N 

All 
All 
0-60W 
60-90W 

41.58(116) 
39.39(87) 

49.87(28) 

11.97 (40) 
12.51(34) 

.43(3) 

*(   )  -  Indicates  sample  size. 
**At  least  source  or  receiver  below  mixed  layer. 


*** 


Both  source  and  receiver  in  the  mixed  layer. 
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TABLE  X. 


Acoustic  Data  Records 


Field 

Columns 

Format 

1 

1-2 

I  2 

2 

3-4 

I  2 

3 

5-6 

I  4 

4 

7-8 

5 

9 

A  2 

6 

10 

7 

11-12 

I  2 

8 

13-14 

I  2 

9 

15-16 

I  2 

10 

17-18 

I  2 

19-21 

I  3 

11 

22 

I  1 

12 

23-25 

I  3 

13 

26-28 

I  3 

14 

29-31 

I  3 

15 

32-34 

I  3 

16 

35-37 

I  3 

17 

38-40 

I  3 

18 

41 

I  1 

19 

42-43 

I  2 

20 

44 

I  1 

21 

45-46 

I  2 

22 

47-48 

I  2 

23 

49-50 

I  2 

24 

51 

I  1 

25 

52-53 

I  2 

26 

54-55 

I  2 

27 

56 

I  1 

28 

57-60 

I  4 

Description 

Cruise  Number 

Station  Number 

Hours  (GCT) 

Minutes  (GCT) 

Sign  of  Time  Zone 

Time  Zone 

Day 

Month 

Year 

Range  (thousands  of  yards) 

Range  (hundreds  of  yards) 

Receiving  Ship 

South  depth  (ft.) 

Receiver  depth  (ft.) 


db) 
db) 
db) 
db) 


Propogation  loss 

Propogation  loss 

Propogation  loss 

Propogation  loss 

Sea  State 

BT  Pattern  Code 

Sea  State 

BT  Pattern  Code 

Degrees  (latitude 

Minutes  (latitude 

N  or  S 

Degrees  (longitude) 

Minutes  (longitude) 

E  or  W 

Water  depth  (fathoms) 


2o2  KC 
8  Kc 
16  Kc 
25  Kc 


Ship  1 
Ship  2 
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TABLE  XI. 
Input  Array  Description  -  NK(I) 

I  Description 

1  Cruise  number,  Zero  implies  all  cruises,  otherwise 
specify  cruise  number  desired. 

2  Station  number.   Zero  specifies  all  stations, 
otherwise  specify  station  number  desired. 

3  Time  (GCT) „   Zero  implies  all  times,  otherwise 
specify  lower  bound  on  time  interval  desired. 

4  Time  (GCT) .   Same  as  NK(3)  except  specify  upper  bound 
on  desired  time  interval. 

5  Time  zone.   Zero  implies  all  time  zones,  otherwise 
specify  desired  time  zone.   If  negative  time  zone, 
fill  in  first  with  negative  sign,  otherwise  leave 
first  character  blank. 

6  Montho   Zero  implies  all  months,  otherwise  fill  in 
lower  bound  no  month  interval  desired. 

7  Month.   Same  as  NK(6)  except  specify  upper  bound  on 
month  interval  desired0 

8  Yearc   Zero  implies  all  years,  otherwise  specify 
year  desired. 

9  Range.   Zero  implies  all  ranges,  otherwise  specify 
lower  bound  on  range  interval. 

10  Range.   Same  as  NK(9)  except  specify  upper  bound. 

11  Receiving  ship.   Zero  implies  both  ships,  1  implies 
Ship  NOo  1,  2  implies  Ship  NoQ  2. 

12  Source  depth.   Specify  source  depth  desired. 

13  Receiver  depth.   Specify  receiver  depth  desired. 

14  Latitude.   Zero  implies  all  latitudes.,   This  is  a 
packed  word  which  expresses  latitude  in  degrees 
and  thousandths  of  degrees.   For  example:   40.817 
degrees  will  be  put  in  the  array  as  40817,  and 
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TABLE  XI.   (Cont'd) 


14  31.7  as  31700.   This  will  always  be  the  lower 
(cont'd)   limit o   If  working  with  a  n/n  or  S/S  situation, 

enter  the  smaller  of  the  two  latitudes.   If 
working  with  a  n/S  situation,  enter  the  north 
latitude. 

15  N  or  S.   This  will  always  be  an  N  for  the  N/S 
situation. 

16  Latitude.   Zero  implies  all  latitudes.   This  is 
a  packed  word  (see  NK(14) ) .   This  will  always  be 
the  upper  bound  on  the  latitude  interval.   If 
working  with  a  N/N  or  S/S  situation,  enter  the 
larger  of  the  two  latitude  limits.   If  a  N/S 
situation,  enter  the  south  latitude. 

17  N  or  So   This  will  always  be  S  for  N/S  situation. 

18  Longitude.   Zero  implies  all  longitudes.   This 
is  a  packed  word  (see  NK(14) ) .   This  will  always 
be  the  lower  limito   If  working  with  E/E  or  W/W 
situation,  enter  the  smaller  of  the  two  longi- 
tudes o   If  E/W  situation,  enter  the  east 
longitude 0 

19  E  or  W.   This  will  always  be  E  for  e/W  situation. 

20  Longitude.   Zero  implies  all  longitudesc   This  is 
a  packed  word  (see  NK(14) ) .   This  is  always  the 
upper  limit  on  longitude.   Enter  the  larger  of 
the  two  longitudes  here  if  a  E/E  or  W/W  situation, 
Enter  west  longitude  for  E/w  situation0 

21  E  or  W0   If  E/W  situation,  enter  W  here. 

22  Water  depth.   -1  implies  all  depths,  otherwise 
specify  lower  limit. 

23  Water  depth.   Same  as  NK(22)  except  specify  upper 
limit. 

24  Transmitter  frequency.   13,  14,  15  and  16  imply 
2.2,  8,  16,  and  25  kc  respectively. 

25  Number  of  stations  in  cruise  5. 

26  Range  differential  between  pairs  of  observations. 
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TABLE  XI.   (Cont'd) 

27  Number  of  stations  cruise  7. 

28  Number  of  stations  cruise  8. 

29  Number  of  stations  cruise  9. 

30  Number  of  stations  cruise  10. 

31  Number  of  stations  cruise  11. 

32  Number  of  stations  cruise  120 

33  Enter  first  cruise  number  to  be  used. 

34  Specify  range  break  point. 

35  Number  of  repetitions  of  program  desired.   Must 
have  a  complete  set  of  input  data  cards  for  each 
iteration. 
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TABLE  XII. 


Description  of  Arrays 

Arrays  Description 

A  Stores  the  time  difference  within  pairs  of 

observations  by  cruises 0 

Al  Accumulates  time  differences  over  all  cruises. 

AT  Stores  time  differences  for  pairs  of  observa- 

tions that  have  both  receiver  and  transmitter 
neither  above  or  below  the  thermocline0 

AU  Stores  time  differences  for  pairs  of  observa- 

tions that  have  both  receiver  and  transmitter 
above  the  thermocline. 

AV  Stores  time  differences  for  pairs  of  observa- 

tions with  both  transmitter  and  receiver  below 
the  the rmoc line o 

B  Stores  the  difference  between  adjusted  trans- 

mission losses  within  pairs  of  observations 
by  cruises. 

Bl  Accumulates  adjusted  transmission  loss 

differences  over  all  cruises. 

BT  Same  as  AT  for  transmission  loss  differences. 

BU  Same  as  AU  for  transmission  loss  differences. 

BV  Same  as  AV  for  transmission  loss  differences. 

DBl         Decibel  loss  for  Case  I. 

* 

DB2         Decibel  loss  for  Case  II. 

DB3         Decibel  loss  for  Case  III. 

* 

DB4  Decibel    loss    for   Case    IV. 

* 

DBO         Combines  db  losses  for  Cases  III  and  IV0 

* 

DB  Combines  db  losses  for  Cases  I  and  II. 
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TABLE  XIIo   (Cont'd) 


IDB 


IK 


IN 


LLX1 


LLX2 


LLX3 


LLX4 


MK 


MMK 


N 

NK 
NY 

PDB 
PRNG 
RNG1 
RNG2 


Stores  first  and  second  digits  for  BT  pattern 
codeD 

Array  MK  is  loaded  into  this  array  one  record 
at  a  time  for  computations  in  subroutine  COMP. 

Loaded  sequentially  with  records  from  the 
acoustic  data  tape  in  the  main  program. 

Establishes  midpoints  of  five-minute  intervals 
in  subroutine  SIGMA. 

Stores  BT  pattern  code  for  Ship  No.  1  and  1st 
data  set  in  pair. 

Stores  BT  pattern  code  for  Ship  No.  2  and  1st 
data  set  in  pair. 

Stores  BT  pattern  code  for  Ship  No.  1  and  2nd 
data  set  in  pair. 

Stores  BT  pattern  code  for  Ship  No.  2  and  2nd 
data  set  in  pair. 

Stores  records  from  acoustic  tape  that  have 
the  same  cruise  no. 

Temporary  storage  for  acoustic  tape  records 
when  a  different  cruise  number  other  than  the 
one  being  considered  is  encountered.,   Will  be 
placed  in  the  MK  array  after  computations  are 
completed  for  the  previous  cruise. 

Counts  entries  in  each  five-minute  interval  in 
subroutine  SIGMA. 

The  input  parameters. 

Counts  entries  in  each  15-minute  interval  in 
subroutine  SIGMA. 

Stores  DBO  for  plotting  at  end  of  each  cruise. 

Stores  RNGO  for  plotting  at  end  of  each  cruise, 

Stores  ranges  for  Case  I. 

* 

Stores  ranges  for  Case  II. 
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TABLE  XII.   (Cont'd) 


RNG3 

RNG4 

RNG 

RNGO 

SLOPE  1 

SLOPE  2 

SLOPE  3 

SLOPE  4 

X 

XX 

XRNG1 

XRNG2 

XRNG3 

XRNG4 

XRNGO 

XRNG 

XY 

XZ 


Stores  ranges  for  Case  III. 
Stores  ranges  for  Case  IV.* 

Stores  ranges  for  Cases  I  and  II. 

* 

Stores  ranges  for  Cases  III  and  IV. 

Stores  ranges  for  Cases  III  and  IV. 

Stores  slopes  for  Case  II  by  station  number. 

Stores  slopes  for  Case  III*  by  station  number. 

Stores  slopes  for  Case  IV  by  station  number. 

Stores  Var(Y(   t) )  within  five-minute  intervals 

Stores  top  and  bottom  of  the rmoc line 
information  for  pairs  of  observations. 

Used  in  plotting  RNGl. 

Used  in  plotting  RNG2. 

Used  in  plotting  RNG30 

Used  in  plotting  RNG4. 

Used  in  plotting  RNGO. 

Used  in  plotting  RNG 

Used  in  computing   Var(Y(At))  in  subroutine 
SIGMA. 

Var(Y(At))   within  each  15-minute  interval 
in  subroutine  SIGMA. 


fRefer  diagram  on  Page  39. 
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TABLE  XIII. 


Description  of  Variables 


Variable  Description 

IA        Set  equal  to  the  station  number.   Used  to 
check  for  the  end  of  a  stationc 

UK       Set  to  the  cruise  number  under  consideration. 
Used  in  finding  LL. 

IQ        Counter  for  arrays  collecting  data  over  each 
cruisec 

JJK       May  be  1  or  zero.   When  set  to  zero  by  the 

program,  it  stops  the  iteration  at  the  end  of 
the  cruise  specified  by  NK(1) . 

JJL       A  switching  mechanism  that  allows  computations 
to  proceed  for  cruise  12  at  the  end  of  the 
acoustic  data  file. 

KNT       Specifies  the  number  of  iterations  within  each 
run  of  the  program. 

KK        Cruise  number. 

L         Sets  the  transmitting  frequency. 

LL        Specifies  the  number  of  stations  in  a 
particular  cruise. 

LLL        Allows  for  zeroing  certain  variables  in  sub- 
routine COMP-,  at  the  beginning  of  each 
iteration. 

. LXl       The  BT  pattern  code  for  Ship  1  and  1st  set 
in  data  pair0 

LX2        The  BT  pattern  code  for  Ship  2  and  1st  set 
in  data  pair. 

LX2        The  BT  pattern  code  for  Ship  2  and  1st  set  in 
data  pair. 

LX3        The  BT  pattern  code  for  Ship  1  and  2nd  set 
in  data  pair. 

71 


TABLE  XIIIo   (Cont'd) 

Variable  Description 

LX4        The  BT  pattern  code  for  Ship  2  and  2nd  set 
in  data  pair. 

MM         Counter  for  the  MK  arrayQ 

MMM        Counter  for  the  MMK  array. 

MX         Counter  for  arrays  collecting  data  over 
all  cruises. 

NC         Set  to  the  cruise  number  under  consideraion. 
Used  to  check  for  the  end  of  a  cruise. 

NN         Sets  the  station  number  for  use  with  the 
SLOPE  arrays. 

SIG  Standard  deviation  of  Y(   t) 

VAR  Var(Y(   t)  ) 

X  Source  depth 

Y  Receiver  depth 

Z  Maximum  allowable  range  differences  between 

pairs  of  observations. 
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